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Bump-on-tail instability of twisted excitations in rotating cold atomic clouds 
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We develop a kinetic theory for twisted density waves (phonons), carrying a finite amount of orbital 
angular momentum, in large magneto optical traps, where the collective processes due to the exchange of 
scattered photons are considered. Explicit expressions for the dispersion relation and for the kinetic (Landau) 
damping are derived and contributions from the orbital angular momentum are discussed. We show that for 
rotating clouds, exhibiting ring-shaped structures, phonons carrying orbital angular momentum can cross 
the instability threshold and grow out of noise, while the usual plane wave solutions are kinetically damped. 
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Introduction. It is an established fact that electromag¬ 
netic radiation can carry a finite amount of angular mo¬ 
mentum, due to the resulting contributions from the spin 
and orbital angular momentum degrees of freedom. Al¬ 
though the spin is always an intrinsic variable, orbital 
angular momentum (OAM) may be either intrinsic or ex¬ 
trinsic [1]. Spin to orbital angular momentum conversion 
may occur in an homogeneous and isotropic medium [2]. 
More recently, attention has been given to the exchange 
of orbital angular momentum between EM radiation and 
different systems, such as electrostatic waves in plasmas 
[3] , nonlinear media [4] , acoustic vortices in optical fibers 
[5] or even sub-meter sized metallic objects [6]. Never¬ 
theless, OAM is not a feature exhibited only by light. 
In fact, vortex phonons have also been proposed to ma¬ 
nipulate particles and small objects [7, 8]. In magnetic 
crystals, for example, spin-orbit coupling can promote 
the excitation of phonons carrying a finite amount of or¬ 
bital angular momentum [9]. 

The manipulation of OAM became attractive in cold 
atomic systems as well since the observation of optical 
pumping of OAM in Cs atoms [10] and the vortex nu- 
cleation in Bose-Einstein condensates [11], by means of 
Laguerre-Gauss (LG) beams. Cold atoms confined in 
large magneto-optical traps are particularly interesting, 
as the atomic dynamics is significantly altered due to 
the appearance of a collective interaction mediated by 
the multiple scattering of light [12, 13]. In this regime, 
the atoms experience a Coulomb-like long-range interac¬ 
tion [14], allowing for an analogy with a one component 
trapped plasma [15]. 

In this paper, we theoretically describe a combined ef¬ 
fect of OAM and the multiple-scattering of light: the on¬ 
set of a dynamical instability in rotating clouds of cold 
atoms triggered by the OAM of the hybrid phonons - 
the elementary excitations in large magneto-optical traps 
(MOT). Hybrid phonons, with properties similar to both 
plasma and acoustic waves [16], have been at the origin 
of a large set of exotic phenomena such as phonon-lasing 



FIG. 1. (color online) Illustration of the experimental setup, 
(a) The transverse beams (darker red), which are slightly mis¬ 
aligned, are of much higher intensity than the beams along the 
longitudinal direction (lighter red), (b) Perspective and (c) 
top views of the density profile for a rotating trap. The cen¬ 
tral cloud (blue) is surrounded by atoms located at the shell 
of mean radius R and width w rotate at the angular velocity 
O, as a consequence of the misalignment of the trapping laser 
along the transverse direction, (d) Poloidal velocity distri¬ 
bution function. The speed distribution of the atoms of the 
rotating shell contains a positive derivative (bump-in-tail), 
which will originate the instability. 


[17]. classical rotons [18], photon bubbles [19] or the dy¬ 
namical Casimir effect [20]. In our setup, and contrary 
to what happens with light - where the phase singular¬ 
ities are imposed by phase modulation - the collective 
excitations (phonons) in cold atoms acquire OAM due to 
the rotation of the cloud, as a consequence of the total 
angular momentum conservation. 

The experimental setup we have in mind consists of 
a large laser-cooled s 'Rb gas, with a typical total num¬ 
ber of atoms N ~ 10 9 — 10 11 . Via a slight misalignment 
of the trapping laser beams along, as described in the 
experiments by Sesko et al. [12], we produce a Saturn- 
like structure (see Fig. 1). The external clump of atoms 
then rotates at a constant angular velocity Ll of the or¬ 
der the thermal velocity of the effective plasma frequency 
oj p [16]. The modification we require in respect to usual 
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experimental configurations is that the laser beam in¬ 
tensity along the longitudinal direction should be much 
smaller than the transverse ones, such that cloud is cigar 
shaped. As a result, the equilibrium of cooling and trap¬ 
ping produce the configuration of a cigar-shaped Saturn- 
like structure as illustrated in Fig. 1. The rotation of the 
exterior shell transfer angular momentum to the excita¬ 
tions (phonon) on top of the equilibrium configuration. 
We therefore show that such “twisted” phonons lead to 
a kinetic instability, a type-1 0 oscillatory instability [21], 
in a rotating cloud of cold atomic trap operating in the 
multiple-scattering regime. This instability is formally 
very close to the bump-on-tail instability observed in flu¬ 
ids and plasmas [22, 23], as the rotation induces an off¬ 
set in the poloidal velocity distribution, allowing for the 
instability to develop. The important difference with re¬ 
spect to other system stems in the fact that only twisted 
phonons can trigger the instability in our case, while zero 
OAM excitations are usually enough for the bump-on-tail 
sort of kinetic instabilities. 

Kinetic description. The complete kinetics of cold 
atoms in magneto-optical traps is described by a Vlasov- 
Poisson equation of the form 

+ v • V + ^ • V„) W (r, V, t) = 1[W( r, v, t)} (1) 

where W (r, v, t) is the distribution function normalized 
to the total number of particles N = f dr f dvW (r, v, t ). 
The total force F = F tra p + F c accounts for the cooling 
and trapping terms, F tra p — —/cr — av, with k and a be¬ 
ing the spring constant and the cooling rate, respectively, 
and the collective force F c = — W c [14-16, 24] with 

V 2 V C (r, t) = —Q J dvW (r, v, t) (2) 


and (2), with W = Wq+SW and V c = V c p+8V c , yielding 

^ + v ■ V ) SW = — ■ V v W 0 , and (3) 

ot ) m 

V 2 SV C = —Q J dv SW (v). (4) 

In what follows, we consider a cigar-shaped trap elon¬ 
gated along the 2 —direction, in the absence of rotation 
(see Fig. 1). Experimentally, this aspect ratio is obtained 
by decreasing the intensity of the trapping lasers along 
the ^-direction. In this case, a slowly-varying amplitude 
density fluctuation propagating along the 2 —axis as e lkz 
must be solution of the following equation 

V 2 SV C ~ - k 2 + 2ifc SV C . (5) 

Under this approximation, the potential SV C also sat¬ 
isfies the paraxial equation (Vj_ + 2 ik-^) 8V = 0 [3]. 
Such a condition reduces Poisson’s Eq. (4) to k 2 8V c = 
Q f dv8W (v). The solution to the paraxial equation can 
be given as a linear combination of Laguerre-Gauss (LG) 
modes as 

SV C (r, t) = J2 Sv pI f pI (A z) e iW e ikz - iut (6) 

P i 

with SVpi being the mode amplitude (assumed to be 
small) and the integers p and l the radial and azimuthal 
mode numbers, respectively. The orthogonality of the 
LG modes reads ( F p i \ F p n /) = 8 pp '8w, with the mode 
functions F p i (r, z) given by 

F pl (r,z)=C pl ^L^e-^ 2 (7) 


where Q = — &l)I/c represents an effective 

charge of the atoms, with ctr and cfl the scattering and 
absorption cross sections [12, 13, 25, 26], respectively, 
and I the light intensity. The term X[W] in Eq. (1) rep¬ 
resents the collision integral and describes the source of 
hot reservoir atoms and the collisions between them [27]. 
As we are interested in small deviations from equilibrium 
only (thus ruling out far-from-equilibrium effects) which 
occur at much slower time-scales that the dynamics of hot 
atoms, we may safely neglect the contribution of the col¬ 
lision integral. Moreover, for large traps (typically with 
N ~ 10 9 — 10 10 ), and appropriate laser detuning, the 
density profiles are approximately constant [28, 29], or 
at least slowly varying at the scale of the fluctuations 
we are about to describe. As such, we should invoke 
the local-density approximation (LDA) which amounts 
to neglect the trapping inhomogeneities. 

In order to account for small amplitude fluctuations 
around the LDA equilibrium Wo, we linearize Eq. (1) 


with £ = r 2 /w 2 (z ) and w(z) the wave beam waist. 
The normalization constants C p i are given by C p i = 


l 


b+p)! 

p! 


1/2 


with the usual Laguerre polynomials, 

4(0 = ^[C‘ + ^]. 

The fluctuations of the collective force are determined 
by (5F C = — ik 0 ffSV c where the effective (paraxial) wave 
vector k 0 ff reads [30] 


k c ff 


i dF p i 
F p i dr 


l 

e r H —eg + 
r 



i dF p i \ 
F p i dz J 


e z . 


( 8 ) 


Similarly, the fluctuations SW can also be decom¬ 
posed into a superposition of LG modes as 8W (v) = 
E pl SW pl (yr)F pl (r,z)e iie e ikz - iut . By plugging in Eq. 
(3), multiplying by F p i and making use of the orthog¬ 
onality relations, we obtain 


6W p i = 


c,pl 

to (a + ib) 


V„Wo, 


(9) 
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with 6F c , p i = -q e ffi 5V p i and q cff = -iq r e r + lq g e e + 
( k — iq z ) e z . Here, we set 

poo p\j-p f. OO 

Qj= / F p i—^rdr, qg = F 2 t dr ( 10 ) 
Jo a J Jo 

with j = {r, z}. The real quantities a and b are defined 
as a = (oo — kv z ) — Iqevg and b = ( q r v r + q z v z ). Notice 
that no mixture between the model p and l is expected 
at the linear level, which is the sufficient description to 
understand the nature of low lying excitations. From 
Eq. (9) and the decomposition of 5W in LG modes, and 
noting that a + ib = (oo — q G ff • v), we finally derive the 
kinetic dispersion relation 


1 - 


oo: 


p 

k 2 


qcff * ^ i; Wo 
(oo - q e ff • v) 


dv=l + x(w) = 0, 


( 11 ) 


where we have defined the quantity Wq(v) = 
Wo(v, r)/no(r) as normalized to the LDA density rio(r), 
and the effective plasma frequency oo 2 = Qn^/m [15, 16]. 
This kinetic dispersion relation encodes the dynamics of 
twisted (l ^ 0) excitations above the equilibrium. De¬ 
pending on the details of the velocity distribution W (v), 
the OAM momentum may change the stability conditions 
of the system, as shown below. 

Twisted Landau damping. Assuming that \q r \, |g~| -C 
\lqg\, and separating the susceptibility into its real and 
imaginary part x( u ) = x!( u ) + ix"(oo), we rewrite Eq. 
(11) in the more explicit way 


( 12 ) 

where V denotes the Cauchy principal part of the inte¬ 
gral. The two poles in the integration occur respectively 
for u z = (to — Iqgvg) /k and ug = (oo — kv z ) /Iqg. On the 
other hand, for the imaginary part of the susceptibility 
we can write 


*' M = | 


dW 0 /dv z 
(V z - u z ) 


dv + V 


8Wo/dvg 
(vg - Ug) 


9Wo\ + fdWp\ 
dv z ) Vz=Uz v dl, o ) v 9 =u 9 

( 13 ) 

For low enough temperatures, it is plausible to assume 
excitations with a phase speed much larger than the 
width of the distribution, u z v z and ug vg, which 
allows us to expand the denominator in Eq. (12) as 


X 


( \ p 

h = -'f 


V 


/ 


dW 0 dvj 1 
7 - t dv ~ — 

{Vj - Uj) u 2 


Wo (Vj) 



dv, (14) 


with j = {z,6}. The real part of the mode frequencies 
can be readily obtained by setting x ' (oo) = 0, which read 


2 2 
u = to p 



3 cik 2 


(15) 



FIG. 2. (color online) Real (left panel) and imaginary (right 
panel) part of the dispersion relation for g 2 = 2 (gray line) and 
g 2 = 5 (red dashed line). The usual phonon dispersion (cor¬ 
responding to the plane wave limit g 2 oo) [15, 16] (black 
solid line) is presented here for comparison. 


where we have defined the speed of sound as c 2 = (i> 2 ) = 
{ V D + ( v g) ir with (vf) = f v 2 W 0 (v) dv and j = {z, 0}. 
The paraxial parameter 77 , accounting for the OAM of 
the excitation, is given by 77 = k/lqg. Assuming that 
axial propagation is dominant, g > 1 - where the paraxial 
approximation makes sense - we get a small correction to 
the hybrid phonon modes derived in [16]. 

Let us now examine the kinetic (Landau) damping of 
the twisted excitations. For that purpose we go back to 
Eq. (13) and note that the damping rate, 7 , is given 
by 7 = — oux" (oo)- The overall damping rate 

is the sum of two contributions, 7 = 7 Z + 7 g, from the 
longitudinal and perpendicular subspaces, respectively. 
We can now assume a Maxwellian distribution of the form 
Wo ,j ( Vj) = l/c s v / 7 re _ D / ' c =, with j = {z,6} and finally 
obtain 

_ _ 2 y/mop 

1 ~~(kX D ) 3 

(16) 

with Ad = c s /oOp being an analogue of the Debye length, 
denoting the typical size of the hybrid phonons in the 
trap. In the limit of zero OAM, g 2 —> 00 , we recover the 
expected result 7 = — e ~ 1 /( fcAp ) 2 . These results, 

as well as the dispersion relation given by Eq. (15) are 
illustrated in Fig. (2). As we can observe, 7 is negative 
for all modes, which corresponds to the Landau damp¬ 
ing mechanism. However, excitations carrying OAM are 
shown to be less damped than the plane wave solutions. 
In what follows we shall demonstrate that, in the case 
of a rotating cloud of cold atoms, and for a particular 
range of achievable experimental conditions, the vortex 
phonons can cross the instability threshold and grow out 
of noise while the plane wave solutions remain kinetically 
damped. 

Bump-on-tail instability. Rotating clouds of cold 
atoms, resulting from small misalignments in the trap¬ 
ping beams, often present ring-shaped density profiles 
[12, 31-33], which are well described by a radial density 
function n (r) = / w , where R is the ring ra¬ 

dius and w its width (see Fig. 1). In this case, the system 


e -l/(k\o) 2 e -l/ri 2 + ve -ri 2 /(^n) 2 e -V 2 
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can be approximately described by a rotating rigid body, 
vg (r) = fir [ 12 ], where fi is the angular frequency of rota¬ 
tion. By plugging this condition in the previous density 
function, we can derive a normalized poloidal velocity 
distribution function 

Wo M = 2 / Ac *, (17) 

Acg^TT 


with Acg = wfi denoting the width of the distribution. 
The interpretation of the later as a distribution func¬ 
tion is justified by the typical scales of the experiments 
[12]. For a rotating frequency of fi ~ (27 t)100 Hz and a 
ring width w ~ 0.5 mm, we get A eg = u>fi ~ 0.3 m/s. 
Such a velocity distribution width is of the same order 
of the thermal velocity of a 81 Rb MOT at T = 200 /xK 
(c s ~ 0.25 m/s). A typical ring radius of R ~ 2 mm cor¬ 
responds to a poloidal velocity of fi R ~ 1 m/s. Plugging 
the new poloidal velocity distribution in Eq. (13), and 
noting that we can approximate A d = c s /lo v ~ A cg/uj p , 
which significantly simplifies the calculation, yields a new 
expression for the kinetic damping term 


7fi =- 


( k\o ) 3 


e -l/(k\o) 2 e -l/v 2 + 


7 ] ( 1 - A (fcA D ) ) e-V/ffcAx,) 2 e -Tr*( 1 -/5/„) a 


,( 18 ) 


with (3 = fi/c s . Note that, for j3 = 0, which corresponds 
to a regular Maxwell-Boltzmann distribution for non ro¬ 
tating clouds, we recover the previous result of Eq. (16), 
as expected, where small oscillations in the homogeneous 
solution are always damped. These results are illustrated 
in Fig. 3 and stability diagrams in terms of the appropri¬ 
ate parameters, q 2 and /3, are sketched in Fig. 4 and Fig. 
5. A type I Q (l/0,w/ 0) [21] oscillatory modulational 
instability can then occur in rotating clouds, where the 
rotation velocity can be understood as an external con¬ 
trol parameter. Experimentally, the rotation velocity can 
be controlled via the misalignment of the trapping beams 
and their intensity. 

In conclusion, we have developed a kinetic theory for 
elementary excitations - hydbrid phonons - carrying or¬ 
bital angular momentum in large magneto optical traps, 
taking into consideration the collective effects emerging 
from the exchange of scattered photons. The resulting 
kinetic dispersion relation, formally similar to that of a 
plane wave, was shown to depend on an effective wave 
vector, gr e ff, in which the helical structure of the wavefront 
is embedded. Explicit expressions for the dispersion re¬ 
lation and kinetic damping rates were derived for a cold 
gas in thermal equilibrium. The contribution from the 
orbital angular momentum states was encoded into the 
dimensionless paraxial parameter q = k/lqg , which is ba¬ 
sically the ratio between the longitudinal and azimuthal 
wave vectors, k, and qg , respectively. By setting q —> oo 
we recovered the expected result for a plane density wave. 




FIG. 3. (color online) Left panel: Kinetic damping rate for 
rj 2 = 2 and different poloidal velocities, /3 = 2 (black line) 
and /3 = 3 (red line). The shaded region corresponds to 
an instability. The gray line corresponds to the non-rotating 
cloud, (3 = 0. (right panel) Kinetic damping rate for constant 
poloidal velocity, /3 = 2 , and different paraxial parameters, 
rj 2 = 2 (red line) and q 2 = 5 (black line). The plane wave 
solution, q 2 —> oo corresponds to the gray line. 



P P 

FIG. 4. (color online) Stability diagram for vortex phonons 
in the (/3, rj 2 ) plane for k\o = 2 (top left), k\D = 2.5 (top 
right), kXo = 5 (bottom left) and kXo = 10 (bottom right). 
The coloured region corresponds to the range of parameters 
for which a modulational instability occurs. Blue regions cor¬ 
respond to higher growth rates. 


The case of a rotating cloud of atoms in a ring-shaped 
structure was considered, in which a new poloidal ve¬ 
locity distribution function was introduced and justified 
based on the typical experimental scales for the system 
temperature, radius and angular velocity. A positive 
damping term from the poloidal velocity distribution can 



FIG. 5. (color online) Stability diagram for vortex phonons in 
the (/3, fcAo) plane for q 2 = 3 (left panel) and in the (t? 2 , fcAo) 
plane for (3 = 2 (right panel). 
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be large enough to compensate for the negative axial 
damping, in realistic experimental scenarios. In this way, 
vortex phonons with a finite orbital angular momentum 
were shown to cross the instability threshold while the 
plane waves remain kinetically damped. This opens the 
possibility of easily detecting this vortex instability in 
rotating clouds of cold atoms in large magneto optical 
traps. Note that the instability mechanism resembles 
the “bump-on-tail” instabilities usually found in plasma 
physics scenarios, with non-equilibrium velocity distri¬ 
butions where regions of positive slope, dWo/dv > 0, 
act as free energy source for the waves, feeding the in¬ 
stability. In our case, the “bumpy-tail” velocity distri¬ 
bution emerges from the rotation of the cloud and, by 
scanning through the f3 parameter, a kinetic instability 
should be triggered and easily detected through stan¬ 
dard fluorescence, or absorption imaging techniques. As 
a final remark we note that the ring-shaped structure is 
not mandatory for the instability to occur, as the essen¬ 
tial ingredient is always a rotation and a “bump” in the 
poloidal velocity distribution, where angular momentum 
from this region is transferred to the twisted wave. In 
order to observe the onset of instability in such a situa¬ 
tion, a filtered image of the fluorescence of the outer shell 
should be recorded. A real-time analysis should therefore 
contain oscillations in the fluorescence signal, at the fre¬ 
quency of the most unstable mode. The measurement of 
the most unstable mode can then be performed by tak¬ 
ing multiple realizations of the Fourier transform of the 
fluorescence signal. 
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